First measurement of coherent 0-meson photoproduction on deuteron at low energies 
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The cross section and decay angular distributions for the coherent <j> meson photoproduction on 
the deuteron have been measured for the first time up to a squared four- momentum transfer t = (p 7 — 
P4 ,) 2 = -2GeV 2 /c 2 , using the CLAS detector at the Thomas Jefferson National Accelerator Facility. 
The cross sections are compared with predictions from a re-scattering model. In a framework of 
vector meson dominance, the data are consistent with the total </>-N cross section ct^jv at about 
10 mb. If vector meson dominance is violated, a larger ct^jv is possible by introducing larger f-slope 
for the <j>N —* 4>N process than that for the ■yN — > (f>N process. The decay angular distributions of 
the <j> are consistent with helicity conservation. 

PACS numbers: 13.25.-k; 13.75.-n; 14.40.Cs ; 25.20.Lj 



The exchange of gluons between hadrons, known as 
Pomeron exchange |l|, is a fundamental process that is 
expected to dominate hadron-hadron total cross sections 
at high energies. In general, multi-gluon exchange is 
harder to study at lower energy since diagrams including 
quark exchange play a more important role. The <f> me- 
son is unique in that it is nearly pure ss and hence multi- 
gluon exchange is expected to dominate </>-N scattering at 
all energies. Since gluon exchange is flavor blind, infor- 
mation on multi-gluon exchange, isolated from the <p-N 
interaction, would be universal and useful in models of 
hadron-hadron interactions. For example, information 
on the <p-N interaction at very low energies, known as 
the QCD van der Waals interaction, is essential for the 
reliable prediction of the possible formation of a bound 
state in the </>-N system [2J. 

The total </)-N cross section (ct^at) is estimated by us- 
ing vector meson dominance (VMD) applied to exclusive 
4> photoproduction on the proton in the photon energy 
range E n < 10 GeV, resulting in o^n— 10-12 mb [EIJ], 
which is in agreement with the estimate from the addi- 
tive quark model [f| applied to KN and ttN scattering 
data [6|. More recently, the inelastic </>-N cross section 
a ™N was extracted from the attenuation of (/>-mesons in 
photoproduction from Li, C, Al, and Cu nuclei 0. The 
attenuation for large A is significantly larger than that 
calculated from VMD. More sophisticated models 0,(1, @| 
are consistent with the experiment if a^jf is significantly 
larger (~30 mb) compared with u^n from the VMD 
model. The reason for the discrepancy of ct^n from these 
two estimates is not well understood. Here we will show 



that information on the t dependence and spin structure 
of the (jyN interaction provides essential clues to solve 
this problem. 

In this Letter, the 4>-N interaction is investigated in 
coherent photoproduction on deuterium. The diagrams 
of the dominant processes contributing to the reaction 
7(f — ► 4>d are shown in Fig. [T] In the first diagram, 
Fig. Q^a), the <f> is produced in a single scattering of a 
nucleon, which is dominant at small —t, and strongly 
suppressed at larger —t due to the deuteron form factor. 
The second diagram, Fig.QJb), shows double-scattering, 
where the <fi is produced at the first vertex and scat- 
ters from the other nucleon at the second vertex. The 
strength of the second interaction is gauged by ct^at. 
The probability to undergo double-scattering increases 
at larger —t because both nucleons receive momentum 
transfer and may recombine into a final-state deuteron 
with a smaller relative momentum between the two nu- 
cleons 0. 

The <f> meson is a spin one particle which decays to a 
KK pair, i.e. two spin-less particles. The decay angu- 
lar distribution of the 4> carries information on the spin 
structure of the reaction amplitude which is the sum of 
single- and double- scattering processes [111 ], 

The measurement of the differential cross sections of 
coherent <j> photoproduction and the decay angular dis- 
tributions in a wide t range allows one to study the <^-N 
interaction in both single and double scattering, as well 
as the transition from one to the other. 

The data were collected with the CLAS detector and 
the Hall B tagged-photon beam at the Thomas Jefferson 
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FIG. 1: (a) Single-scattering and (b) double-scattering con- 
tributions to the coherent </>-meson photoproduction on the 
deuteron. 



National Accelerator Facility [l2j. The incident electron 
beam energy was 3.8 GeV, producing tagged photons in 
the range from 0.8 to 3.6 GeV. The photon beam was 
directed onto a 24-cm long liquid-deuterium target. The 
data acquisition trigger required two charged particles 
detected in coincidence with a tagged photon. Charged 
particles were momentum analyzed by the CLAS torus 
magnet and three sets of drift chambers. The torus mag- 
net was run at two settings, low field (2250 A) and high 
field (3375 A), each for about half of the run period. 

The reaction 7c? — > <f>d was identified by detecting 
a deuteron and a K + from <p K + K~ decay. The 
K + and deuteron were selected based on time-of- flight, 
path length, and momentum measurements. Figure [2]^ a) 
shows the missing mass distribution, Mx, for the reac- 
tion 'yd — ► dK + X when events near the 0-meson peak 
(0.98 < M(K+K-) < 1.12 GeV/c 2 ) were selected in the 
K + K~ invariant mass, assuming a K~ was the miss- 
ing particle. A missing K~ peak is seen on top of a 
smooth background from non-<i K + K~ final states. The 
missing mass resolution, ranging from 8 to 30 MeV/c 2 , 
depends on photon energy and the deuteron momentum. 
A three-cr cut was applied to select the missing K~ for 
the exclusive jd — > K + K~d reaction. 

Figure [2fb) shows the invariant mass distribution 
for the K + K~ pair after the selection of the missing 
K~ . The 0-meson peak appears above a smooth back- 
ground. The 0-meson yield was obtained from a fit 
to the M(K+K~) distribution by a gaussian-convoluted 
Breit-Wigner function and a background function. The 
width and the pole position for the Breit-Wigner func- 
tion were fixed to 4.3 MeV/c 2 and I0I9.5 MeV/c 2 , re- 
spectively [lij]. The standard deviation of the gaus- 
sian distribution was fixed to the value obtained from 
simulation. The background function was chosen as 
b{x 2 



a^Jx 2 — (2m K ) 2 + b(x 2 — (2rriK) 2 ) [lJ], where x is 
M(K + K~), mx is the charged kaon mass, and a and 
b are the fit parameters. Three background functions: 
a linear background, background from non-resonant 
K + K~d production, and /q photoproduction, were stud- 
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FIG. 2: (a) Missing mass distribution of the reaction 
— > dK + X for events near the 0-meson mass (0.98 < 
M(K + K~) < 1.12 GeV/c 2 ). (b) Invariant mass distribution 
for the K + K~ pair after the selection of the missing K . The 
solid curve is a fit to the data. The dashed curve shows the 
contribution from background. 



ied as alternative choices. The background models for the 
non-resonant K + K~d and /o photoproduction were pa- 
rameterized by the differential cross section and photon- 
energy distribution of events in the sidebands of the 4>- 
meson peak. The dependence of the yield on the back- 
ground function, fit range, and parameterization of the 
Breit-Wigner function were studied. The extracted yield 
changes between 3% and 9% depending on the yield ex- 
traction procedures. 

The CLAS acceptance was determined by using a 
GEANT-based Monte Carlo simulation [l5j]. The simula- 
tion was iterated to reproduce measured t, photon energy, 
and decay angular distributions. The acceptance was be- 
tween 10% and 20% in the kinematic region covered by 
the present measurements. The accuracy of the calcu- 
lation of the acceptance was estimated from comparison 
of results from the other event reconstruction topologies 
(d K + K~, K + K~ . and d topologies) for which the 
acceptances were different from that for the d K + topol- 
ogy. The differential cross sections for these topologies 
are shown in Fig. [3J They agree with each other within 
statistical uncertainties, indicating that the acceptance 
is understood as to the number of reconstructed tracks, 
charge combinations, and decay modes. Supplemental 
simulations were performed to understand the system- 
atic uncertainties due to the event generator (1-11%) and 
event reconstruction (1-5%). 
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TABLE I: Differential cross sections for the reaction 7c? — > <j>d. 
The second and third numbers in each field are the statistical 
and systematic uncertainties, respectively. 
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FIG. 3: Comparison of differential cross sections for 'yd — > <f>d 
from various topologies in the range 1.6 < E- t < 3.6 GeV. 
Only statistical uncertainties are shown. 



Systematic uncertainties in the yield extraction and ac- 
ceptance were estimated as a function of photon energy 
and t; they were between 4% and 13%. The combined 
systematic uncertainty for the luminosity and trigger ef- 
ficiency was less than 10%. Systematic uncertainties 
from contributions from accidental tracks, target win- 
dows, and particle misidentification are less than a few 
percent. The total systematic uncertainty was estimated 
as 11-17% by adding these uncertainties in quadrature. 

The differential cross sections were measured in the 
ranges 1.6 < E 1 < 2.6 GeV and 2.6 < E 1 < 3.6 GeV. 
They are given in Table U Figured] shows the experimen- 
tal data in the range 2.6 < E 1 < 3.6 GeV. The data is 
compared with theoretical calculations using a rescatter- 

13, 



ing model [IOJ, [16[ . In this model, the 7./V — > <j>N ampli- 
tude was parameterized by using published data on the 
7P —> 4>p reaction [171 ] and data from the proton target 
run during this experiment. This amplitude was convo- 
luted with the deuteron wave function with a correction 
of the relativistic-recoil effect [lOj . The double scattering 
process (Fig.QJbVF is modeled by the Generalized Eikonal 
Approximation [18]. The er^jv and t dependence for the 
re-scattering process are the inputs for the calculation. 
The model successfully reproduces the differential cross 
sections on coherent p photoproduction [Io| using the in- 
puts from the VMD. 

The total model uncertainty is estimated to be 
about 20%. A 10% uncertainty was assigned to the 
parametrization of the 7-/V — > <f)N amplitude based on 



the 7p — > 4>p data. The effect of spin-flip in the process 
7 A" — > <pN was ignored in the parametrization of the 
single scattering amplitude since the spin-flip amplitude 
is more suppressed in the coherent process than in the 
incoherent process. A 15% systematic uncertainty was 
assigned due to this effect [19(. An isospin dependence 
of the process 7^ — > <j> N was not taken into account in 
the model, but Ref. [20J suggests such an effect is small. 

In Fig. |H curve A shows the t distribution calculated 
by using the VMD prediction for the 0-N cross section, 
i.e. er0Ar=lO mb, and the same t distribution for the reac- 
tion 7 A" — > (f)N and the reaction tpN — > cf>N. The curve B 
corresponds to a c f > N=30 mb, inspired by Ref. with the 
VMD assumption for the t distribution. It overestimates 
the data at large —t where the contribution from double 
scattering dominates. This implies that if the t distri- 
bution follows the VMD prediction, ct^n should also be 
consistent with the VMD prediction. In this case, in- 
consistency with the larger a^N from the ^-dependence 
experiment 0] still remains. 

However, the VMD picture may not be a good ap- 
proximation in this photon energy range. The larger 
(70 n from the ^4-dependence experiment [3] can be ex- 
plained if the t distribution of the reaction — » cf)N 
is different from the VMD prediction. For example, it is 
possible for the virtual <f> to fluctuate to a KK pair and 
have a larger cross section for the second interaction 21 1 . 
In this case, the i-slope for the second interaction would 
be larger than that for the 7^ — * <j)N reaction based on 
a general geometric relation between the t-slope and the 
total cross section (22[. Following this hypothesis, cross 
sections were calculated with cT0jv=3O mb using a larger 
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FIG. 4: Differential cross sections for the reaction "fd — » 4>d. 
The inner error bars shown are statistical uncertainty only, 
while the outer error bars are the sum of statistical and sys- 
tematic uncertainties in quadrature. The curves A, B and C 
are calculations from the re-scattering model [Tol . , see text 
for details. The uncertainties on curves B and C are compa- 
rable to that of curve A, but are not shown. The dot-dashed 
curve is a contribution from the single scattering diagram. 



exponential i-slope, b^N = 10 (GeV/c)~ 2 , in the second 
interaction (curve C). The data is equally- well described 
by curve C, suggesting a larger t-slope parameter is neces- 
sary if (JcfiN is larger than the VMD prediction. Although 
the current data do not allow one to extract the o^n and 
the i-slope independently due to the strong correlation 
between them, it possibly suggests a larger ct^n from the 
A-dependence. 

In addition to the differential cross sections, the decay 
angular distributions of the <f>. meson were also measured 
in the helicity frame The direction of the 4>- meson 
momentum in the CM system was chosen as the z-axis, 
and the polar angle and azimuthal angle between the 
K + momentum and the 0-meson production plane were 
defined as Oh and 4>h in the </>-meson rest frame. Fig- 
ure [5] shows the projections of the decay angular distri- 
butions onto cos6h and 4>h in the ranges —0.8 < t < 
-0.35 GeV 2 /c 2 and -2.0 < t < -0.8 GeV 2 /c 2 in each 
photon energy region. The data are consistent with the 
prediction from helicity conservation (solid curves), i.e. 
the spin of the </>-meson is aligned to the momentum of 
the (i-meson. This is similar to what was observed in the 



smaller —t to the larger —t region, implying that the spin 
structure of the single- and double-scattering processes 
are similar. 

In summary, we have presented the first measurement 
of the differential cross sections and decay angular distri- 
butions for coherent <fi photoproduction on the deuteron 
up to t = —2.0 GeV 2 /c 2 . The differential cross sections 
at large —t exhibit a contribution from double scattering. 
The data are consistent with er0jv=lO mb in a framework 
of VMD. The data also provide a possible explanation for 
larger ct^n if the t-slope for <fiN —> <pN is larger than the 
VMD value from yp —t (fyp. The decay angular distribu- 
tions follow the prediction from helicity conservation. 

This measurement demonstrates a new approach to 
the study of the <fi-N interaction in the low energy re- 
gion where VMD is not necessarily a good approximation. 
Further measurements at higher photon energies [25| . at 
very small —t [2a|. as well as an A-dependence study in 
e + e~ decay [271 ] will make it possible to map out details 
of the energy and t dependences of the <^-N interaction. 
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(f) photoproduction on the proton [23j, [24( . In the larger 
—t region, the double scattering contribution becomes 
more important. No drastic change is observed from the 



tt 
[1] 

m 

[3] 



Current address:University of New Hampshire, Durham, 
New Hampshire 03824-3568 

Current address:TRIUMF, 4004, Wesbrook Mall, Van- 
couver, BC, V6T 2A3, Canada 

Current address:Massachusetts Institute of Technology, 
Cambridge, Massachusetts 02139-4307 
Current address:Catholic University of America, Wash- 
ington, D.C. 20064 

Current address:Physikalisches Institut der Universitaet 
Giessen, 35392 Giessen, Germany 

Current address:Edinburgh University, Edinburgh EH9 
3JZ, United Kingdom 

Current addressUniversity of Massachusetts, Amherst, 
Massachusetts 01003 

A. Donnachie, H. G. Dosch, P. V. Landshoff, and 
O. Nachtmann, Pomeron Physics and QCD (Cambridge 
University Press, 2002). 

H. Gao, T. S. H. Lee, and V. Marinov, Phys. Rev. C63, 
022201(R) (2001). 

H. J. Behrend et al., Nucl. Phys. B144, 22 (1978). 



6 
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